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Abstract 

We derive the lepton spectrum in semileptonic beauty decays from a nonperturbative treat- 
ment of QCD; it is based on an expansion in 'i/mq with mg being the heavy flavour quark 
mass. The leading corrections arising on the l/^i^Q level are completely expressed in terms 
of the difference in the mass of the heavy hadron and the quark. Nontrivial effects appear in 
l/mg terms affecting mainly the endpoint region; they are different for meson and baryon 
decays as well as for beauty and charm decays. 

The weak decays of hadrons contain a wealth of information on the fundamental 
forces of nature. Yet the intervention of the strong interactions has prevented us from 
extracting this information in a reliable way. Heavy flavour decays promise to be more 
tractable since the mass of the heavy flavour quark uiq provides a powerful expansion 
parameter. Indeed significant successes have been scored by Effective Heavy Quark 
Theory (EHQT)|l|]. Yet that approach has some intrinsic limitations: e.g., it deals 
with exclusive semileptonic modes only and it requires the presence of heavy quarks 
both in the initial and in the final state; thus it cannot be applied directly io h ^ u 
transitions. Its model-independent predictions so far include corrections through 
order l/mg only. On the other hand the energy released in heavy flavour decays is 
much larger than ordinary hadronic energies. Our analysis will make use of this large 
energy release in treating Q qlv transitions with mg, ttlq — uiq ^ Aqcd- 

In previous papers we have shown how nonperturbative contributions to global 
quantities like lifetimes and semileptonic branching ratios can be obtained from a 
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model-independent treatment of QCD. The method was based on expanding the 
weak transition operator into a series of local operators of increasing dimension with 
coefficients that contain increasing powers of l/mq. The coefficients depend on the 
(inclusive) final state; the differences between the decay rates of different heavy flavour 
hadrons Hq - charged vs. neutral mesons vs. baryons - enter through the matrix 
elements of the local operators taken between Hq. For example the total semileptonic 
h ^ u width through order 1 / is given by: 

V{H, luX) oc ml ■ — ^ {Hb\ bb - \biaGb\H,) (1) 

IMHf, mjj 



with iaG = i'j/j.'juG^u , Gfj,i, being the gluonic fleld strength tensor. 

In this note we will expand the general method to treat the lepton spectra in the 
semileptonic decays of heavy flavour hadrons. A novel feature is encountered: one 
is dealing with an expansion in powers of 1/{pq — PiY with pq and pi denoting the 
momenta of Q and the lepton /, respectively. This series is singular at the endpoint of 
the lepton spectrum; thus some care has to be applied in interpreting the results. Two 
dimension flve operators generate the leading nonperturbative corrections of order 
l/mg : the colour magnetic operator QiaGQ and the operator Q {D"^ — {vD)'^)Q 
describing the kinetic energy of Q in the gluon background fleld; denotes here 
the covariant derivative and is the 4-velocity vector of the hadron. Corrections 
actually arise already on the l/mq level; it is crucial that in QCD those can be 
expressed completely in terms of the difference between the quark and the hadron 
mass. 

We will phrase our discussion in terms of beauty decays with a few added com- 
ments on charm decays. It should be noted that the question of the inclusive lepton 
spectra in QCD was flrst addressed explicitely in ref. 0. 

Ignoring gluon bremsstrahlung one obtains a lenghty expression for the lepton 
spectra in the semileptonic decays of a beauty hadron Hh. 

^{H, IvX,) = T,e{l-y- p) 2y\{l - ff{l + 2/)(2 - y) + {I - ff{l - y) + 

^ in P 

+ (1 - /) [ (1 - /)(2 + -y - 2/ + -fy) - ^{2y + /(12 - 12y + 5y')) ] G, - 

6 6 p 

- [ ^(1 - + 2f)y + - /)(10y - 8y') + ^(3 - 4/)(2y2 - y') ] K, } 
3 p p2 

1927r'^ 1 — y m^ m?, 

K, = l/2MH,-mb{iD)%\H,)lml , ; G^, = II2Mh,- {H^\biaGb\H^) I2ml {2) 

with mq denoting the mass of the quark q in the flnal state. Eq.(2) represents the 
master formula containing both relevant cases, namely q = u, c. 

For b —>■ ulv transitions with m„ = this expression simplifles considerably: 

^ = ro2y2[3-2y-(^y+-<5(l-y)--(2y2-y3)5'(l-y))ir,+ (2+^y--<5(l-y))G,] 
ay 3 3 6 3 6 

(3) 



2 



The (5-functions and their derivative reflect the previously mentioned singular nature 
of the expansion at the endpoint; their emergence has a transparent meaning (see^ 
for details). The spectrum is finite at the endpoint for = and thus contains a 
step function. The chromomagnetic interaction effectively 'shifts' the spectrum by 
changing the energy in either initial or final state; the shift in the argument of the 
step function thus yields a 5-function. The singular structure in the Ki, term on the 
other hand reflects the motion of the b quark inside the Hj, hadron which Doppler 
shifts the spectrum; in second order it generates 6'{1 — y) for the step-like spectrum. 

Due to these singular terms the expression given above can be identified with the 
observable spectrum only outside a finite neighbourhood of the endpoint region. (This 
distance remains constant in absolute units in the limit mt — * oo). Yet even this 
neighbourhood does not represent true 'terra incognita': for integrating our expres- 
sion over this kinematical region yields a finite and trustworthy result that can be 
confronted with the data. This can be expressed through the function 

fEmax 

T{Ei) = dEi — , Ei< Eqcd < E^ax • (4) 

J El aJ^i 

Emax denotes the maximal kinematically allowed energy and Eqcd the maximal en- 
ergy for which one can still trust the QCD expansion in eqs.(2,3); its value depends 
on the size of Kh and Gf,. Clearly r(0) = Tsl has to hold. This is not a trivial 
relation: for r5j;^ is deduced from a completely regular expansion in l/rUb (see eq.(l)) 
whereas r(0) comes from integrating the expression in eq.(3) containing singularities; 
thus the singular terms are essential for recovering the correct decay width! 

Similar considerations apply to 6 ^ clu transitions. With rUc as an infrared cutoff 
there arise no singular terms at the endpoint y = 1 — p; yet the expansion parameters 
Gb/ (l — y), Kh/{1 — yY though finite become large there. The need for 'smearing' the 
spectrum over the endpoint region, eq.(4), thus still exists. 

The contributions from higher dimensional operators that we are ignoring have 
terms of the schematical form ~ [/^hadr/("^6(l ~ y))T] summing them all up would 
yield a well-behaved function. As long as these quantities are smaller than unity, we 
can trust the expressions given above for the unintegrated lepton spectrum. 

The size of G^ is easily determined: G^ = {B\Q iaGQ \B) /imf ~ 3{M'^{B*) - 
M^{B))/Aml ^ 0.017; for A;, it vanishes. A recent QCD sum rules analysis yields |^ 
{B\b {iDf h\B) ~ 0.4 GeV^ ■ 2Mb in agreement with rather general expectations. For 
our subsequent discussion we will set = 0.02. We use m?, ~ 4.8 GeV as deduced 
from a QCD analysis of the Ypsilon system, mb — nic = 3.35 GeV as inferred from the 
beauty and charm meson masses which for mi, = 4.8 GeV yields nic = 1.45 GeV, and 
put m„ = 0. For b ^ u decays we estimate that Eqcd ~ 0.9 ■ mb/2 ^ 2.15 GeV; for 
b ^ c a. somewhat smaller 'smearing' range near the endpoint seems to be required, 
namely AEi ~ 0.15 GeV and thus Eqcd — 2.0 GeV. Numerically it implies that in 
the real world the c quark is relatively heavy in b decays: > /Uj^adr^^b ! therefore 
the falling edge of the spectrum starts in the calculable region. 

For a proper perspective we show the partially integrated spectrum T{Ei) from 
three different prescriptions, namely (a) our QCD expansion; (b) the simple free 
quark picture with Gh = = 0, m„ = 0; (c) the phenomenological treatment of 
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Altarelli et al. (hereafter referred to as ACM) where one attempts to incorporate 
some bound state effects. We have set here m„ = nispect = 0.15 GeV, nic = 1.67 GeV, 
Pf = 0.3 GeV as suggested by a fit to CLEO data; pp denotes the 'Fermi momentum'. 
To be consistent we have ignored gluon bremsstrahlung both in ACM and in our QCD 
expressions (for fitting the data it should be added to eqs.(2-4)). In comparing the 
QCD formula with the ACM prediction one has to note a subtle distinction in the 
definition of the kinematical variables: in ACM energy is expressed in terms of the 
mass of the beauty hadron; yet in eqs.(2-3) y measures the energy in units of mi,, 
which by itself introduces a shift of order l/m^. 

Comparing the results on T{Ei) from the three approaches we conclude: (i) The 
shape of the QCD, of the free quark model and of the ACM curves are very close 
over most of the range of Ei < Eqcd for the b —>■ u as well as the b —>■ c case. 
The main difference lies in the overall normalization, (ii) The QCD result for b —>■ 
u can be largely simulated by setting m„ ~ 0.3 GeT^ in the free quark spectrum. 
The nonperturbative corrections thus effectively transform a current quark mass into 
a 'constituent' one of reasonable size, (iii) Once the ACM result is renormalized 
according to r(0)QCD(^ ^ u) = 1.07 x r(0)ACM(^ ^ u) and r(0)QCD(& ^ c) = 
1.13 X r(0)ACM(^ ~^ c) the two curves are hardly distinguishable as functions of Ei. 
(iv) While the QCD results for T{Ei) are largely independent of the value of K}, for 
Kb ~ 0.02 in 6 — i> M, they are sensitive to it for b ^ c; that change can be easily 
understood|Q. (v) The QCD curves for semileptonic A^ decays are somewhat harder 
than for B meson decays, (vi) The largest differences between the models are found 
in the endpoint region, namely for Ei > 1.8 GeV in 6 — > c. 

Some of these points are illustrated in the figures. In fig.l we show T{Ei) for the 
endpoint region of b ^ u; in fig.2 we have plotted these partially integrated spectra 
for b ^ c. The unintegrated spectrum dT{Ei) / dEi which can be calculated in the 
QCD treatment for Ei < 2 GeV is shown in fig. 3. The similarity in the shape of the 
two spectra is quite striking! 

A few comments are in order about charm decays. The lepton spectra are quite 
different in beauty and in charm decays, especially in the endpoint region: for the 
charged lepton in beauty decays is an electron or muon whereas it is an antifermion 
Z"*" in charm decays. Since the lepton spectrum vanishes at the end point in c — > qlu 
even for = 0, the chromomagnetic operator can yield finite terms only while the 
kinetic energy operator produces a 5(1 — y), but not a 5'(1 — y) term. Yet since the 
nonperturbative corrections are much larger in charm than in beauty - Gc, ~ 0.2 

- the smearing required in y is much larger for charm, namely Ay ~ 0.5. While this 
allows us meaningful statements about fully integrated quantities hke FgL, it appears 
beyond the scope of our present analysis to treat lepton spectra in charm decays. 

To SUTriTnarize: We have shown here how an expansion in l/{pb—pi) allows us to 
incorporate nonperturbative corrections to the lepton spectra in inclusive semileptonic 
b c and b ^ u decays. Due to the singular nature of this expansion at the endpoint 
the spectrum cannot be computed in detail in a close neighbourhood of the endpoint, 
yet smeared or partially integrated spectra can. Through order 1/m^ these expressions 
are given in terms of the quantities m^, m^, Gb and Kb. These are not free parameters 

- their size can be extracted from the relationship to other observables. Improvements 
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in the precision of their determination, in particular concerning Kb and to a lesser 
degree m^, can be anticipated from future progress in theory. 

We have compared our QCD results with the ACM description that so far has 
allowed a decent fit to the data. Despite the obvious differences in the underlying 
dynamics we have found the shape of the resulting lepton spectra remarkably similar, 
and even more so when one keeps the following in mind: 

(a) Initially the scale for the kinematics is different in the two descriptions: for the 
QCD expansion it is set by the quark mass whereas for ACM by the hadron mass 
Mb or M\^. The kinematical differences due to 7^ Mh^ are actually of order 1/mh 
and formally represent the leading corrections. This suggests an interesting observa- 
tion: an accurate measurement of the shape of the spectrum allows a determination 
of the b quark mass rrib free of theoretical uncertainties. Unfortunately in practice one 
would have to analyse b ^ u; the shape of the b —>■ c spectrum is basically determined 
by mfc — rric with little sensitivity to m^. 

(b) The ACM ansatz contains three free parameters - rric, mgp and pp - that are 
to be fitted from the data. There are four quantities that set the scale in the QCD 
expansion through order 1/ml, namely nib, Gb and Kb but none of them is a free 
parameter. That our QCD description containing therefore very little 'wiggle room' 
can nevertheless yield a description so similar to that of the ACM ansatz - after the 
latter has been fitted to the data - has to be seen as quite remarkable! 

(c) The ACM prescription can thus be viewed as a simple, though smart approxima- 
tion to a more complete and complex QCD treatment. At the same time it would be 
incorrect to interprete the fit parameters in ACM literally as real physical quantities; 
it is thus not surprising that the numerical values for the former differ significantly 
from what is known now about the corresponding quantities in QCD. 

A more detailed analysis shows 0] that QCD provides a natural 'home' for the 
Fermi motion originally introduced phenomenologically in ACM, and it is asymptot- 
ically a dominant effect for the end point shape in b decays. However it enters in a 
somewhat different form. Its impact on the total widths is only quadratic in 1/mb 
(in ACM it produces a linear shift). Nevertheless the shape of the lepton spectrum 
gets l/nib corrections^; it will be discussed in detail elsewhere 

Despite all similarities there arise relevant numerical differences as well, not only 
in the overall normalization of the spectra, but also in the endpoint region: the 
QCD spectra place a higher fraction of 6 — >■ c events in the endpoint region. Such a 
difference has important consequences for the extraction of iK^b/Kbl from the data. 
Furthermore our analysis shows that the lepton spectra in semileptonic A;, decays are 
distinct from those in B decays. These issues will be discussed in a future paper [Q. 

Another model ISGW[^ used for describing semileptonic spectra relies heavily on 
model-dependent calculations of the formfactors for the exclusive final states. In the 
limit of a heavy c quark, — -C rric , EHFT yields the necessary formfactors to 
some accuracy; therefore in that limit the prediction of this model coincides with the 
model-independent QCD spectrum. On the other hand that is apparently not the 
case for 6 — >■ m decays where the difference is significant. 

On the l/rriQ level that was not treated here novel effects arise due to 'Weak 

^At this point we disagree with the conclusion of ref.|Bl 
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Annihilation' in the h ^ u channel; for details see 
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Figure Captions 

Fig.l: The function T{Ei) for 6 — > m transitions calcul ated in the QCD expansion 
(sohd hne), the free spectator quark model (dotted hnc) and the ACM ansatz (dashed 
hne); the ACM curve has been multiphed by a factor 1.07. 

Fig. 2: The function T{Ei) for 6 — > c transitions calcul ated in the QCD expansion 
(solid line) , the free spectator quark model (dotted hne) and the ACM ansatz (dashed 
line); the ACM curve has been multiphed by a factor 1.13. 

Fig. 3: The lepton spectrum dV/dEi for b ^ c transitions calculated in the QCD 
expansion (solid line), the free spectator quark model (dotted line) and the ACM 
ansatz (dashed line); the ACM curve has been multiplied by a factor 1.13. 
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